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Abstract 

We present a high-precision QED calculation, with 0.1% theoretical accuracy, of 
two photon production in e + e~ annihilation, as required by more and more ac- 
curate luminosity monitoring at flavour factories. The accuracy of the approach, 
which is based on the matching of exact next-to-leading order corrections with a 
QED Parton Shower algorithm, is demonstrated through a detailed analysis of the 
impact of the various sources of radiative corrections to the experimentally relevant 
observables. The calculation is implemented in the latest version of the event gen- 
erator BabaYaga, available for precision simulations of photon pair production at 
e + e~ colliders of moderately high energies. 
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1 Introduction 



The precision measurement of the hadron production cross section in e + e~ 
annihilation at flavour factories, such as $, r-charm and -B-factories, requires a 
detailed knowledge of the collider luminosity [l| . It can be derived by counting 
the number of events of a given reference process and normalizing this number 
to the corresponding theoretical cross section Q. It follows that, in order to 
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maintain small the total luminosity error given by the sum in quadrature of the 
relative experimental and theoretical uncertainty, the reference process must 
be a reaction with high statistics and calculable with an accuracy as high 
as possible. For this reason, the luminosity monitoring processes employed at 
flavour factories are QED processes, namely Bhabha scattering, two photon 
and muon pair production. In particular, at DA$NE, VEPP-2M and PEP-II 
the large-angle Bhabha process is primarily used and the other reactions are 
measured as cross checks 0,0,0], while at CESR [5] all the three processes 
are considered and the luminosity is derived as an appropriate average of the 
measurements of the three QED reactions. 



At DA$NE, a comparison between the luminosity measurement using Bhabha 
events and the process e + e" — > 77 shows very good agreement, the average 
difference in a run-by-run comparison being 0.3% |6|]. This precision necessarly 
demands progress on the theory side, since the Monte Carlo (MC) programs 
used for the simulation of photon pair production, i.e. BabaYaga v3.5 [7|j8| and 
BKQED [9], have a theoretical precision of about 1%. Actually, the original 
formulation of BabaYaga is based on a QED Parton Shower (PS) approach 
for the treatment of leading logarithmic (LL) QED corrections and, as such, it 
lacks the effect of O (a) non-log contributions, which are important to achieve 
a precision at the per mille level. On the other hand, the generator BKQED 
relies on an exact O (a) diagrammatic calculation, therefore neglecting the 
contribution of higher-order LL corrections, which have been already demon- 
strated to be necessary for O (0.1%) luminosity monitoring at flavour factories 
f3l.f7l.riol] . Because of this motivation, the aim of the present paper is to describe 
a high-precision calculation of photon pair production in QED, based on the 
matching of exact next-to-leading-order (NLO) corrections with the QED PS 
algorithm, along the lines of the approach already developed for the Bhabha 
process in Ref . 10] . This will allow a reduction of the theoretical error in lumi- 
nosity measurements at flavour factories, as demanded, in addition to precision 
measurements of the hadronic cross section, by improved experimental deter- 
minations of the e + e~ — > r + r~ cross section at low energies [ll||, important 



~e~ — > t + t~ cross section at low energies 
for precision calculations of the anomalous magnetic moment of the muon. 
Furthermore, a precise knowledge of e + e~ annihilation in two photons is of in- 
terest for estimates of the background to neutral meson production. We do not 
include in our calculation pure weak corrections, which have been computed 
in Ref. 12j and turn out to be important at very high energies, well above the 
energy range explored by flavour factories. For completeness, it is worth men- 
tioning that an independent calculation, including exact O (a) contributions 
supplemented with higher-order LL terms through collinear QED Structure 
Functions 13j, of the relevant corrections to e + e _ — ► 77 at moderately high 
energies was performed in Ref. 141 ]. recently revisited in Ref. 
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The outline of the paper is as follows. In Section 2 we describe the matching 
algorithm for the e + e _ — > 77 process, while in Section 3 we provide numeri- 



2 



cal results, both for integrated cross sections and differential distributions of 
experimental interest, in order to discuss the effects of the various sources of 
radiative corrections and provide evidence for the per mille accuracy of the 
approach. Conclusions and possible perspectives are drawn in Section 4. 



2 Theoretical formulation 

The cross section of the photon pair production process, with the additional 
emission of an arbitrary number of photons, can be written in the LL approx- 
imation as follows 

OO 1 

d< L = n 2 (Q 2 ,e) £ - |^f n , LL | 2 d$ n , (1) 

n=0 n - 

where II (Q 2 , e) is the Sudakov form factor accounting for the soft-photon (up 
to an energy equal to e in units of the incoming fermion energy E) and virtual 
emission, e is an infrared separator dividing soft and hard radiation and Q 2 is 
related to the energy scale of the hard-scattering process. In our calculation, 
Q 2 is fixed to be equal to the squared centre of mass (cm.) energy s, by 
comparing with the exact O (a) calculation of Ref. ||. \^ n ,Lh\ is the squared 
amplitude in LL approximation describing the process with the emission of n 
additional hard photons, with energy larger than e in units of E, with respect 
to the lowest-order approximation e + e~ — > 77. d$ n is the exact phase space 
element of the process (divided by the incoming flux factor), with the emission 
of n additional photons with respect to the Born-like final state configuration. 
The Sudakov form factor, which is defined as 

n(Q 2 ,e) = exp , (2) 

where 

L = log % J+ = dz V (z) , (3) 
m Jo 

appears in Eq. ([1]) to the second power to account for the presence of two 
charged particles in the initial state. In Eq. ([3]) V(z) is the electron — > electron 
+ photon splitting function V{z) = (1 + z 2 )/(l — z). 

The cross section as calculated in Eq. ([1]) has the advantage that the photonic 
corrections, in LL approximation, are resummed up to all orders in perturba- 
tion theory. On the other hand, the weak point of Eq. ([1]) is that its O (a) 
expansion does not coincide with the exact O [a) (NLO) result. Actually, we 
have 
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a n- 1 Q 



2 ' 



<Kl = I 1 - -2+ In ^ ] |^o| 2 d$ + |^i, LL | 2 d$i 

= (l + C QiLL )|^o| 2 d$ + |^ 1 , LL | 2 d$ 1 , (4) 

whereas an exact NLO can be always cast in the form 

d(7 a = (l + C aiS v)Ko| 2 d$ + Ki| 2 d$i , (5) 

where the coefficient C Qi sv is equal to the exact squared amplitude of the 
annihilation process, in the presence of soft and virtual radiative corrections 0. 



14j . in units of the exact Born squared amplitude |^#o| , and is the 



exact squared matrix element of the radiative process e + e~ — > 777 [16]. The 
matching of the LL and NLO calculation can be obtained considering the 
correction factors (free of infrared and collinear singularities) 

Fsv = 1 + (C a , S y - Cq,ll) F h = 1 + I^! 2 -|^l| 2 (6) 

As can be seen, the exact O (a) cross section as in Eq. ([5]) can be expressed, 
up to terms of O (a 2 ) and in terms of its LL approximation, as 

da a = F sv (1 + C a , LL ) |^o| 2 d$ + F H \^ ll \ 2 d^ . (7) 

A similar procedure, repeated to all orders in a, leads to the correction of 
Eq. (JTJ), which becomes 

00 1 / n \ 

d< atched = ^s V n 2 (Q 2 , e) E ry II *h,< I^,ll| 2 d$ n , (8) 

n=0 \i=0 / 

where 



^ = 1+ '^' 1 , 2 i>LLl , (9) 

with |^#j| 2 and I-^ll) 2 squared matrix elements, exact and in the LL approx- 
imation, respectively, relative to the emission of the i-th hard bremsstrahlung 
photon. The expansion at O (a) of Eq. (jSJ) coincides now with the exact NLO 
cross section of Eq. (jSJ) and higher-order LL contributions are the same as in 
Eq. p. 



3 Numerical results 



3.1 Integrated cross sections: technical tests and radiative corrections 



The calculation of QED corrections requires the introduction of the unphys- 
ical soft-hard separator e. Therefore, the indipendence of the predictions for 
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Fig. 1. O (a) QED corrected cross section of Eq. (5) as a function of the infrared 
regulator e. The error bars correspond to la Monte Carlo statistics. 

the QED corrected cross section from variation of such a parameter has to 
be proved, for sufficiently small e values. This is succesfully demonstrated, at 
a precision level of ~ 0.01%, in Fig. [1] and Fig. [21 which show the cross sec- 
tion of the photon pair production process, obtained according to the exact 
O (a) cross section of Eq. (jSJ) (Fig. [T|) and to the matched cross section of 
Eq. (jHJ) (Fig. |2j), as a function of e from 10 -2 to 10 -6 . The numerical results 
shown in Fig. [T] and Fig. [2] correspond to the following experimental set up, 
which models, up to a good accuracy, the selection criteria adopted by KLOE 
Collaboration at DA$NE E3l 
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(10) 



where E™ m 



is the minimum energy threshold for the detection of at least two 
photons, ■$ niin > max are the angular acceptance cuts and £ max is the maximum 
acollinearity between the most energetic and next-to-most energetic photon. 



As a further test of the approach, we checked that our results for the NLO 
corrections agree at the 0.1% level with those quoted in Ref. [sj for the exact 
O (a) relative corrections to the totally inclusive e + e~ — > 77 cross section, as 
a function of different cm. energies. 

To quantify the overall impact of QED radiation and, in particular, to evaluate 
the size of QED contributions at different perturbative and precision levels, 
we show in Tab. CD the Born cross section a, the O (a) PS and exact cross 
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Fig. 2. The same as Fig. 1 for the matched QED corrected cross section of Eq. (8). 
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3 


10 
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137.53 


15.281 


1.3753 


(T PS 

u a 


128.55 


14.111 


1.2529 


NLO 

u a 


129.45 


14.211 


1.2620 


a PS 

w exp 


128.92 


14.169 


1.2597 


"exp 


129.77 


14.263 


1.2685 


5 a 


-5.87 


-7.00 


-8.24 


5oo 


-5.65 


-6.66 


-7.77 


<^exp 


0.24 


0.37 


0.51 


iNLL 
u a 


0.70 


0.71 


0.73 


6 NLL 


0.66 


0.66 


0.69 



Table 1 

Photon pair production cross sections (in nb) to different accuracy levels and relative 
corrections (in per cent) for the set up specified in the text. 

section, <r PS and o"a LO , respectively, as well as the all-order PS cross section 
of Eq. (CD) (T^jjJZI and the matched cross section of Eq. (JSj) o"ex P - The results 
have been obtained with experimental selection criteria similar to those of 
Eqs. ( TTUj) . but with cm. energies y/s = 1,3,10 GeV and with the condition 
E™ n = 0.3-^/i. From these cross section values, it is possible to calculate the 
relative effect of various QED contributions, namely the contribution of exact 

(a) radiative corrections, of higher-order corrections in the O (a) matched 

1 The results denoted as <J^L agree with the predictions of BabaYaga v3.5 within 
0.1% accuracy, as we checked explicitly. 
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PS scheme, of exponentiation with respect to the exact O (a) cross section and, 
finally, of non-logarithmic terms entering the O (a) cross section and present 
in the improved PS algorithm. The above per cent corrections are shown in 
Tab. [T]and they can be derived from the cross section values according to the 
following formulae 

s a = 100 x °L Soo = 100 x ffcxp 

a a 

_ NLO NLO _ PS 

f) cxp = 100 X NLO °<x = 100 X PS I 11 ) 

® a. ® a. 

PS 

5 NLL = 10Q x ggg ~ ^cxp 
"exp 

The numerical errors coming from the MC integration are not shown in Tab. [1] 
because they are beyond the quoted digits. From Tab. [Q it can be seen that 
the exact O (a) corrections, measured by the relative contribution 5 a , lower 
the Born cross section of about 5.9% ($ resonance), 7.0% (J/ip resonance) 
and 8.2% (T resonance). All-order corrections, due to the presence of an ar- 
bitrary number of photons and measured by the relative contribution 8^, 
amount to about 5.7% ($ resonance), 6.7% (J/ip resonance) and 7.8% (T res- 
onance), showing that the introduction of higher photon multiplicity gives an 
increasing of the O [a) corrected cross section. Such an effect, due to O (a n L n ) 
(with n > 2) terms, is quantified by the contribution S exp , which is a positive 
correction of about 0.2% ($ resonance), 0.4% (J/i/j resonance) and 0.5% (T 
resonance), and, therefore, important in the light of the aimed per mille accu- 
racy. On the other hand, also next-to-leading O (a) corrections, quantified by 
the contribution 5^ LL , are necessary at the precision level of 0.1%, since their 
contribution is of about 0.7%, almost indipendently of the cm. energy. Their 
effect is unaltered at the level of 0.1% by the matching procedure with PS, 
as can be inferred by comparing <5^ LL with 5^ L . To further corroborate the 
precision reached in the cross section calculation, we also evaluated the effect 
due to the most important sub-leading O (a 2 ) photonic corrections and given 
by a 2 L contributions enhanced by infrared logarithms. Actually, the bulk of 
such corrections is effectively incorporated in our approach, by means of fac- 
torization of O (a) next-to-leading terms with the leading O (a) contributions 



taken into account in the PS scheme, as argued and demonstrated in Ref. [18 
It turns out that the effect due to O (a 2 L) corrections, which can be inferred 
from the cross section values according to the following formula 



a - a NLO _ PS , PS 
r mn CX P Ce u exp ' u a 
0a 2 L = 100 X 



does not exceed the 0.05% level. 
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As a whole, these results demonstrate that both next-to-leading O (a) and 
multiple photon corrections are unavoidable for 0.1% theoretical precision. 

3.2 Differential distributions 

In Fig. [3] and Fig. H] we show the angular and energy distribution of the most 
energetic photon, while in Fig. [5] the acollinearity distribution of the two most 
energetic photons is represented. The above distributions, which have been 
simulated by using the latest version of the generator BabaYaga, correspond 
to the experimental set up of Eq. ffTUj) and refer to exact O (a) corrections 
matched with the PS algorithm as in Eq. (jSJ) (solid line), to the exact NLO 
calculation as in Eq. ([5]) (dashed line) and to all-order pure PS predictions of 
BabaYaga v3.5 [s| (dash-dotted line). In the inset of each plot, the relative 
effect due to multiple photon contributions (5 exp ) and non-logarithmic terms 
entering the improved PS algorithm (<5^, LL ) is also shown, according to the 
definitions given in Eq. (TTTT) . 




20 40 60 80 100 120 140 160 

■d (deg) 



Fig. 3. Angular distribution of the most energetic photon according to the PS 
matched with O (a) corrections (Eq. (jSJ, solid line), the exact O (a) calculation 
(Eq. ([5|), dashed line) and the pure all-order PS as in BabaYaga v3.5 (dash-dotted 
line). Inset: relative effect (in per cent) of multiple photon corrections (solid line) 
and of non-log contributions of the matched PS algorithm (dashed line). 

For the angular distribution of the most energetic photon, the contribution 
of higher-order corrections beyond O (a) amounts to about 1% in the cen- 
tral region and is at a few per mille level at the edges of the distribution. 
More pronounced effects due to exponentiation are present for the the energy 
distribution of the most energetic photon. In the statistically dominant re- 
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Fig. 4. The same as Fig. 3 for the energy distribution of the most energetic photon. 

gion around 0.5 GeV, higher-order corrections reduce the O (a) distribution 
of about 20%, while they give rise to a significant hard tail in the proximity of 
the energy threshold of 0.3^, as a consequence of the higher photon multi- 
plicity of the resummed calculation with respect to the fixed O (a) prediction. 
Concerning the acollinearity distribution, the contribution of higher-order cor- 
rections is positive and of about 10% in correspondence of quasi back-to-back 
photon events, whereas it is negative and decreasing from ~ —30% to ~ —10% 
for increasing acollinearity values. 

1000 , , 




Fig. 5. The same as Fig. 3 for the acollinearity distribution of the two most energetic 
photons. 

As far as the contributions of non-logaritmic effects, dominated by next-to- 
leading O (a) corrections, are concerned, they contribute at the level of some 
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per mille for the angular and acollinearity distribution, while they lye in the 
some per cent range for the energy distribution. 



Therefore, as for the cross section, the interplay between O (a) corrections 
and exponentiation is crucial for precise predictions at the level of differential 
cross sections. 



4 Conclusions 



We have presented a high-precision QED calculation of the e + e~ — ► 77 pro- 
cess, of interest for luminosity monitoring at flavour factories. The calculation, 
which includes all the relevant radiative corrections at the 0.1% precision level, 
is implemented in an improved version of the event generator BabaYaga[l], 
in order to contribute to a reduction of the uncertainty in luminosity mea- 
surements at e + e~ colliders of moderately high energies. The accuracy of the 
approach has been demonstrated through a careful analysis of the various 
sources of photonic corrections to the experimentally relevant observables. 
The per mille precision reached for the two photon production process can, 
indeed, take advantage of the absence in such a process of the vacuum polar- 
ization uncertainty, present in both e + e~ — > and e + e~ — > e + e~, and can 
also rely on previous estimates of sub-leading two-loop corrections to Bhabha 
scattering, that have been shown to typically contribute at the level of a few 



0.01% in the energy range explored at flavour factories (see e.g. Ref. jlOj and 
references therein). 

Possible perspectives concern the application of the approach, here presented 
for e + e~ — > 77 and in Ref. [l(| for the Bhabha process, to obtain precise 
predictions for e + e~ — ► and e + e~ — > fi + fi~j, both of interest for physics 

studies at flavour factories, as well as the inclusion of the exact O (a) weak 
corrections to the two photon production process. 
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